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Introduction 

In the past three decades the thermal noise spectroscopy was recognized as a fruitful tool of 
space plasma diagnostics. It is well-known that when a passive electric antenna is immersed into a 
stable plasma, the thermal motion of the ambient electrons and ions produces fluctuations of the 
electric potential at the antenna terminals [Rostoker, 1961; Andronov ,1966; De Passiz, 1969; 
Fejer and Kan, 1969]. This quasi-thermal noise may be calculated if the particle velocity 
distribution function is known [Rostoker, 1961]. Since the noise spectrum depends on the main 
characteristics of plasma, as electron density or bulk temperature, the quasi-thermal noise 
spectroscopy can be used for diagnostics of plasma parameters. This diagnostic method is most 
appropriate for in situ space plasma measurements because it does not require of additional 
apparatus; large scale space plasma conditions allow one to construct antennas whose impedance is 
very small as compared with input impedance and whose characteristics can be calculated with fair 
accuracy. Such antennas permit direct observation of the frequency spectrum of thermal noise. 
Some examples of application of this diagnostic method can be found, for example, in [Meyer- 
Vernet, 1979; Couturier et al, 1981; Kellog, 1981; Meyer-Vemet and Pershe, 1989]. 

It is also well known that random irregularities of electron density always present in real 
space plasma. Random irregularities of the Earth's ionosphere are studied intensively and main 
properties of their spatial spectrum are known [Fejer and Kelley, 1980; Szuszczewicz, 1986]. These 
irregularities considerably affect propagation of radio waves in space plasma changing their phase, 
amplitude, spatial and angular distribution [Zabotin, Bronin and Zhbankov, 1998], as well as group 
propagation time and pulse duration [Bronin, Zabotin and Kovalenko, 1999]. Some information 
about the irregularity spectra of solar wind is also available [Rickett, 1973]. For the purposes of 
present investigation it is possible to use twin models of shape of spatial spectrum for both 
ionospheric and solar wind irregularities. Though that does not relate to parameters of the spectra. 

Irregularities substantially change properties of the medium with relation to 
electromagnetic radiation and they may also influence quasi-thermal noise spectrum detected by 
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antenna. What is the possible physical mechanism of this influence? It is known, that fluctuations 
in plasma are closely connected with plasma dissipative properties. From the viewpoint of 
statistical mechanics, random irregularities in plasma may be understood as non-thermal large scale 
(in comparison with characteristic time and scale of particle motion) fluctuations. Such fluctuations 
may considerably change collision term in kinetic equation [Klimontovich, 1982] and, 
consequently, velocity distribution function. Since the quasi-thermal noise spectrum is determined 
by velocity distribution function, the change in distribution function will lead to change in noise 
spectrum. From the viewpoint of electrodynamics, random irregularities change mean dielectric 
properties of the medium (see, for example, [Hagfors,1984]). It means that in the media with 
random irregularities the roots of dispersion equation are shifted with relation to undisturbed values 
and even new roots may appear. The imaginary part of these roots may also be considerably 
changed due to additional non-coUisional attenuation caused by scattering of waves in random 
medium. 

Development of strict theory of influence of random irregularities on quasi-thermal noise 
spectrum is a very complicated and difficult task, involving various fields of electrodynamics and 
statistical physics. In the present paper we will assume that noise spectrum in the random medium 
may be calculated using the tensor of effective dielectric permittivity. This tensor is determined as 
the dielectric permittivity tensor of some imaginary "effective" regular medium in which the field 
of point source is the same as the mean field in the corresponding random medium [Ryzhov, 
Tamoikin and Tatar skii, 1965, Ryzhov 1967, Ryzhov 1968]. Since effective dielectric permittivity 
tensor determines the mean field of the source in random media it also determines the impedance 
of antenna, which in its turn, determines the noise spectrum. It has been shown also that correlation 
function of electric field fluctuations in random medium may be expressed through the effective 
dielectric constant tensor [Ryzhov, 1968]. (It should be noted, however, that this result is based on 
averaging of Kallen-Welton formula and is valid only for states near the thermal equilibrium.) 
Spectrum of electrostatic noise in its turn is completely determined by correlation function of 
longitudinal electric field [Meyer-Vernet and Per she, 1989]. The method based on effective 
dielectric constant tensor was used to study thermal noise spectrum in dielectrics [Ryzhov, 1970]. 
We apply this approach to the irregular space plasma. 

Calculation of the effective dielectric permittivity tensor of plasma is also rather complex 
problem involving summation of infinite series of perturbation theory. Some approximation is 
necessary to get the tensor components in closed form. In this paper we follow Born 
approximation. 

Thus, in the present paper we intend to estimate possible changes in the thermal noise 
frequency spectrum in plasma with random irregularities rather then to develop complete and strict 
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theory of this effect. In Section 1 of the paper we briefly discuss calculation of the noise spectrum 
and impedance of antenna in Maxwellian plasma (we assume that presence of irregularities does 
not change velocity distribution function radically and it may be approximately described by 
Maxwell distribution). In Section 2 the method of calculation of effective dielectric permittivity 
tensor is stated. Section 3 discusses specific for given problem difficulties in numerical 
calculations. Results of numerical calculations in application to ionospheric plasma and plasma of 
solar wind are collected in Section 4. 

1. Quasi-thermal noise spectrum and antenna impedance 

Usually the noise voltage spectral density measured at the electric antenna terminals, which 
is immersed in a plasma may be expressed through the spatial Fourier transform of the current 
distribution in the antenna and autocorrelation function of electrostatic field fluctuations in the 
antenna frame [Rostoker, 1961]. The shot noise, drift of the plasma across the antenna as well as 
some other phenomena also contribute into noise spectra. We do not take these phenomena into 
account now. If the plasma is in the thermal equilibrium at temperature T , what will be assumed in 
further consideration, the voltage spectral density may be expressed through antenna impedance by 
the formula [Meyer-Vemet and Pershe, 1989]: 



where kg is Boltzmann constant, Z is antenna impedance and Re denotes the real part. 

Calculation of antenna impedance in general case is rather complex problem. One must take 
into account various phenomena such as the disturbance of trajectories of particles, collection of 
electrons and ions, photoemission and so on. The common approximation is to take into account 
only electron plasma oscillations. In this approximation one has 



where E||(a),^) is the plasma longitudinal permittivity, J(k) is the spatial Fourier transform of 

the antenna current distribution. For Maxwellian coUisionless plasma longitudinal dielectric 
permittivity is given by the well-known expression [Akhiezer, 1972]: 
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where 4>(z) = 2zexp(-z^) J"exp(?^)<i? , z = cak ^ J [-sflQi) pk) , kD=l/rj), rj^ is Debye's length, 



(Hp is plasma frequency. 

To calculate impedance one needs to choose the current distribution J(k) . The frequently 

used model of current distribution corresponds to wire dipole antenna [Kuhl 1966; Kuhl 1967; 
Kellog 1981]. The geometry of wire dipole antenna is shown in Fig. 1. The antenna consists of two 
cylinders, each of length L and radius a« L. This model may be a good approximation for 
geometry of real antennas used for space missions. The antenna parameters used in our calculations 
are in close correspondence with real parameters of spacecraft antennae used in investigations of 
ionosphere or solar wind (Wind and Kosmos-1809 missions). 

For wire dipole antenna one can get [Meyer-Vemet and Pershe, 1989; Kuhl, 1966, Kellog, 

1981]: 

4 Imeii(a),fc) 
Re(Z) = — \dk -F(k ) , (3) 
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where 

F(k) = Uo (ka)f ^ ^ , 
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Jq(x) is Bessel function of the first kind and Si(x) is integral sinus function (see, for example, 

[Abramowitz and Stegun, 1964]. To calculate quasi-thermal electromagnetic noise spectrum in the 
plasma with random irregularities one should substitute the longitudinal dielectric permittivity for 

regular plasma e||(a),^) by effective longitudinal dielectric permittivity e^-^(a),^) [Ryzhov, 
1968]. 

2. Calculation of effective dielectric permittivity tensor 

General scheme of calculation of effective dielectric permittivity tensor starts from the 
Dyson equation for the mean field in random medium (see [Rytov, Kravtsov and Tatarskii, 1978]). 
It may be easily found that effective dielectric permittivity tensor is proportional to the Fourier 
transform of the mass operator {Ryzhov, 1967]. Taking into account only initial term in the mass 
operator series expansion (what corresponds to Bom approximation for the mean field) one should 
obtain [Ryzhov, 1967] 



ef ((o,k) = 8°. (to, £) - k^Bi^ij((o, k) Jg^; (CO, p)^ ^(k-p)d^p, (4) 

where B^^y = (e^ -5(„, )(8^- -5y ), 8^ (C0,£) is dielectric permittivity tensor of regular plasma, 
G^i (co,p) is the Fourier transform of Green tensor, kQ =(xi/c, c is velocity of light in vacuum, 

(^) is spatial spectrum of irregularities. 
In isotropic plasma one has 
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where longitudinal permittivity 8||((0,A;) is defined by formula (3) and transversal permittivity is 
determined by the expression [Akhiezer et al, 1972]: 

8j.(to,fc) = l-^(4>(z)-mi/^exp(-z2)). 

At this point of our consideration we should choose a concrete type of the irregularity 

spatial spectrum ^j^(k). Ionospheric irregularities are described by complex spectrum, which 

properties are different for different intervals of wave numbers of irregularities, but for the certain 
interval of wave numbers the spectrum may be described by power law: 

Generally the spectrum is anisotropic: irregularities may be strongly stretched along the 
lines of force of geomagnetic field due to difference in the diffusion coefficients for longitudinal 
and transversal directions. However, for simplicity of calculations, we will use the model of 
isotropic spectrum. This situation is less characteristic of topside ionospheric plasma, but may be 
considered normal for solar wind plasma. 

In isotropic medium with isotropic irregularities tensor of effective dielectric permittivity 
has the same structure as tensor in regular media (5) (in other words effective medium is also 
isotropic): 
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For specific spectrum index |x = 4 (quite typical value both for ionospheric and for solar 
wind plasma) the spectrum can be written as follows: 
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where / = /27t, is the outer scale of spectrum, 5/j = D(R) , D(R) is structure function of 

irregularities at scale length R . Substituting (10) into (8) one obtains for longitudinal effective 
dielectric permittivity: 
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Internal integration in (1 1) can be done in anal5^ical form: 
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These functions turn into zero when ^ is zero (what corresponds to limiting cases 
p— >0,/7— >oo,/c— >0 or A;— >oo) and provide proper convergence of integrals in (11) at both limits 
of integration over p . 

Though one could use result (11) directly, without substitution of analytic expressions for 
the inner integral, generally such approach may be not successful. The reason is that integrands in 
(11) have several peculiarities, which will be discussed below, as well as significant dependence of 
integrand on such parameters as Debye's radius, outer scale of the irregularity spectrum, 
dimensions of antenna etc. The latter difficulty, however, may be considerably reduced by 
introduction of dimensionless variables. 

3. Details of numerical calculations 

Considerable difficulties under numerical calculation of quasi-thermal noise spectrum and 
of effective dielectric permittivity tensor using expressions (3) and (11) correspondingly are caused 
by the "bad" behavior of integrand at some special points. In both cases integrand contains 
expressions of one of the following kinds: 

Im£||(co,^) ^^^^ 
{Ree^^((xi,k)f +{lm^^^((xi,k)f 
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As an example, in Fig. 2 the dependence of Re en (solid lines) and Im8|| (dashed lines) on 
krD for four values of dimensionless frequency CO/ co^ (= 1.01; 1.05; 1.1; 1.2) has been plotted. 

Functions (12) and (14) have a peak of height ~ l/ImEn and width ~ Im8|| (or ~ ijlmef and 

~ Im8|^-^) when Re8|| =0 or Re8|^-^ =0. For co~ co^ the values of Im8|| are very small when 
Re8|| = and the peak is very sharp (see Figs. 4 and 5). Function (13) has two peaks of opposite 
signs at the left and right of the point where Re8|| =0, and these peaks are also very sharp at the 
frequencies close to plasma frequency (see Fig. 3). A usual way to calculate integrals with such 



quasi-singular integrand is to split the integration interval at the point where the singularity 
happens. In opposite case the result would depend upon the position of the peak with relation to 
the grid of abscissas of integration rule which would be different for different frequencies. In other 
words, in this case the calculated noise spectra would contain random error component due to 
accumulation of inaccuracy under numerical integration. However the splitting of integration 
interval is not sufficient for successful integration. Additional difficulties are caused by the fact that 

the roots of equation Re 8|| =0 may be found only approximately ( Re 8|| -10"^^ in the root), while 

the magnitude of ImEn in the root may be several orders of magnitude smaller. To avoid this 

interference one should take into account the fact that in real plasma some small coUisional 
attenuation is always present. The value of collision frequency may be chosen so that inequality 
Re 8 II « Im 8|| always holds in the small vicinity of the root of equation Re 8|| = . 

As it has been noticed above it is useful to introduce dimensionless variables to reduce the 
dependence of numerical calculations on certain values of plasma parameters. It is natural to 

determine dimensionless frequency ro and wave number k as 

k = kQk, co=co^ro. (13) 

In these variables longitudinal dielectric permittivity, for example, may be written as: 
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In (14) plasma properties are taken into account only through the 
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dimensionless constant ^Qr^ . 

Using (12) - (14), expression (11) for effective dielectric constant tensor can be 
transformed to the following form 
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It is easy to show, that the second integral in (15) is at least (kQl) times smaller than first and 
since may be neglected. 



Integrand of first integral in (15) contains factor 



l+{kolf{p-kf 



which also has 



singular-like behavior when kf^l » 1 . This problem may be solved by appropriate substitution of 

independent variable. For example, the substitution p = k +tan'&/kQl transform integrand into 

smooth function without singularities at p = k . 

We used Brent algorithm for searching the roots of equation Re8|| (S, A: ) = ; the adaptive 

integration method based on Gauss-Kronrod rule was used to perform numeric integration in 
expression (15). 

4. Results of numerical calculations 

For calculations of quasi-thermal noise spectrum we first chose the regular parameters of 
plasma and spectrum of irregularities corresponding to the ionosphere F-region: plasma frequency 
(x)p /In was equal to 3 MHz, Debye's length - 5 cm; R = I km, / =10/271 km. Calculations 

were carried out for four different values of5/j: 510 ,10 ,210 and 3 10 . The first 
value is typical for undisturbed mid-latitude ionosphere, the last is observed in polar and equatorial 

ionosphere and in the experiments of plasma heating by the powerful HF radio wave. 

In regular plasma the spectrum has a peak just above the plasma frequency; the shift from 

the plasma frequency is of order of 5o)/ O)~(ro/L)^ [Meyer-Vemet and Per she, 1989]. For the 
chosen parameters of antenna ro/L«land rp/a«l.It means that noise spectrum does not 
depend upon the radius of antenna wire and have a sharp peak near the plasma frequency. The 



dimensionless spectra — V for regular plasma and for irregular plasma with given are 

16^57 



plotted in Fig. 6. 

In presence of irregularities we observe noticeable change in the spectrum of quasi-thermal 
noise. The peak of the spectrum in this case is split into two peaks. The splitting of the peak is 
caused by the complication of dielectric properties of plasma in the presence of random 
irregularities. At the relatively high magnitudes of 5/^ real part of effective longitudinal dielectric 
constant has additional roots, while the contribution of the roots existing in regular plasma is 
dumped due to increase of the imaginary part caused by scattering. The distance between peaks is 
greater for greater irregularity level. Its magnitude is small, approximately 10 - 15 KHz for 
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5jj = 5-10 and -100 KHz for 5^ =3-10 . Such details of the spectrum, however, can be 
detected in experiment. 

For the solar wind plasma we chose the following parameters: plasma frequency 

(Hp 11% = 20 KHz, Vj) = 10 km, / = 10^ / 2% km. Noticeable influence of random irregularities on 

the electromagnetic noise spectrum takes place at relatively high level of irregularities. In Fig. 7 the 

noise spectrum is plotted for 5^ = 5 • 10"^ , 15 • 10"^ and 30 • 10"^ for normalization scale length 

i? = 10^m. The dimensionless noise spectrum for 5jj =15 10"^, 30 10"^ and 50-10"^ for 

normalization scale 7? = 1 m is plotted in Fig. 8. In the case of solar wind plasma we observe 
the formation of plateau in the spectrum in the vicinity of plasma frequency instead of the splitting 
of the peak. This difference is explained by the difference in ratio of Debye's length to wire length 
for ionospheric plasma and for solar wind plasma. 

Conclusions 

In this paper we have considered the influence of random irregularities of electron density 
in isotropic plasma on the quasi-thermal noise spectrum using fairly simple model of irregularities. 
We have found that for the small values of irregularity level modification of noise spectrum is 
negligibly small. However, for larger values of AA^ / N , also quite possible in natural conditions, 
irregularities cause some noticeable effects. In the ionospheric plasma it is the splitting of the peak 
in the frequency noise spectrum located just above the plasma frequency, into two peaks. Though 
the gap between those peaks is small it still may be detected in experiment. The magnitude of the 
gap depends upon the value of AA/^ IN , what makes possible using measurements of noise spectra 
for the purpose of the irregularity diagnostics. In the solar wind plasma irregularities cause changes 
of the shape of the main spectrum maximum near the plasma frequency resulting in appearance of 
the plateau under higher irregularity level. This effect also can provide essential information about 
the solar wind irregularities. 

In this paper we used simplified model of plasma which may be significantly improved. 
The major possible improvement concerns the spectrum of irregularities. For example, in real 
topside ionosphere irregularities are stretched along the lines of force of geomagnetic field and 
their spectrum is anisotropic. One can also take into account drift of the plasma across the antenna, 
because noise spectra measurement are done on the board of satellite moving through the plasma. 
Though the account for magnetic field is considered to be unimportant for calculation of noise 
spectrum in regular plasma, in the presence of irregularities its influence on wave propagation may 
be important. 
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For other kinds of space plasma, like the plasma of solar wind, the account of all these 
factors may be essential, first of all because such plasma may be anisotropic even in the absence of 
external magnetic field. Besides generally it cannot be considered as being in the thermal 
equilibrium, so more general approach to deriving expressions for noise spectra may be required. 
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Fig. 1 . Wire dipole antenna geometry 
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Fig. 2. Real and imaginary components of longitudinal dielectric permittivity of a regular plasma 
at the frequency near plasma frequency. Figures at the curves indicates the ratio of wave frequency 
to plasma frequency, corresponding to the curve. 
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Fig. 6. Dimensionless quasi-thennal noise spectrum in the ionosphere. 
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Fig. 7. Dimensionless quasi-thermal noise spectrum in the solar wind; normalization scale 

R = ^O^m. 
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Fig. 8. Dimensionless quasi-thermal noise spectrum in the solar wind; normalization scale 

R = ^O^m. 



